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ABSTRACT 


Due to the lack of longitudinal stiffening along 
submarine hulls, they are inherently vulnerable to 
underwater explosions directly off the bow or stern. 
Accordion-like deformations of the hull are set into motion 
which could cause dynamic amplification in the transient 
response of attached substructures. In underwater shock 
acceptance tests of internal equipment, this interaction is 
created by exploding a charge in a fore and aft configura- 
tion with the submerged shock test vehicle (SSTV). With the 
inereasime avallability of large computers and the rapid 
development of numerical methods, several computer codes 
have been written to predict equipment response to 
underwater shocks. Using the ELSHOK (ELASTIC SHOCK) code, 
this investigation studies the effect of hull/substructure 
interaction on stiffened shell response at resonance 
following an end-on load. The transient response of the- 
coupled shell/substructure system from tapered and 
conventional charges of equivalent impulse is examined in 


this study 
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lon TRODUC TION 


A. BACKGROUND 

Submarines are inherently vulnerable to explosions 
detonated directly off the bow or stern. While they have 
adequate protection against side loads with transverse 
bulkheads and frames, there is little longitudinal 
protection along the hull. As a result of this, they are 
susceptible to accordion-like deformations of the hull upon 
impact from end-on loads. It 1S possible that the frequency 
of the hull motions could excite internal equipment into 
resonance thereby causing unacceptable damage. This 
investigation examines the coupled elastic response of a 
shell/substructure system when subjected to an end-on 
underwater explosion (UNDEX) as the mass and stiffness of 
the internal substructure is varied. | 

The current specifications for conducting underwater 
shock tests on submarine-installed equipment are contained 
in MIL-S-901D [Ref. 1]. This document specifies the 
explosive charge weight and geometry of the test, as well as 
the mounting and orientation of the equipment being tested 
within the submerged shock test vehicle (SSTV). The SSTV is 
nothing more than a ring-stiffened cylindrical shell with 


Circular endplates designed to simulate the hull motions 


2 


created by an underwater shock. Recognizing the potential 
for unacceptable damage to hull mounted equipment from an 
end-on load, the specification requires the first of four 
UNDEX's to be performed in a fore and aft configuration. 
The suitability of an equipment design or installation is 
evaluated according to its ability to function as intended 
during and after each shock impulse. Equipment tests of 
this nature are quite expensive and require a great deal of 
preparation. It is desirable for the designer to have some 
idea of what the transient response of the hull and 
equipment will be prior to the actual UNDEX. 

The development of analytical methods and computer codes 
for the analysis of the reaction of submerged structures to 
underwater explosions has enabled the engineer to simulate 
the shock response of submarines and equipment with 
increasing accuracy. Pe eeelemeniaiinite difference 
methods allow for the analysis of Structure and-fluid 
responses, and the work of Geers [Ref. 2] provides a means 
PomIlncomeeratomLluta-structure interaction effects. Several 
ComplreecruNecodesscurrently exist which utilize these 
principles to analyze installed equipment response to shock 
waves. One of them is the ELSHOK (ELASTIC SHOCK) code 
developed by Weidlinger Associates, under the sponsorship of 
the Defense Nuclear Agency (DNA) and the Office of Naval 


Research (ONR). This code was developed to investigate 


r3 


modern submarine underwater explosive shock response in 
conjunction with a testing program using small to large 
scale models and shaped or tapered explosive charges. The 
accuracy of ELSHOK has been validated in several highly 
controlled tests, notably the 1983 low level explosive test 
of an SSN 668 class submarine. The code was made available 
to the Naval Postgraduate School by the DNA with support 
from Weidlinger Associates. It was first used at the school 
by LT Mark Welch, USN, in his study comparing the shock 
response predictions of ELSHOK with those obtained using the 


Dynamic Design Analysis Method (DDAM) [Ref. 3]. 


B. PURPOSE FOR THIS INVESTIGATION 

An explosive detonated underwater exerts a great amount 
of pressure on the surrounding fluid. The water is actually 
compressed by the force of the explosion causing a shock 
wave to form at the point of detonation amen propagates in 
a roughly spherical shape at the speed of sound in water. 
When the pressure pulse impinges on a submarine, dynamic 
responses result in the hulle fromeenewt Wand—s eaueeure 
interaction. The submarine experiences translational motion 
away from the point of impact. The velocity of this motion 
is of particular interest to a designer of submarines or 
weapons. 

In an explosion with a side-on aspect with the 


Submarine, the shock wave hits the hull in the middle before 
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it hits the ends due to its spherical shape. This causes 
the middle to bow away from the point of impact while the 
ends remain initially fixed. The resulting motion is called 
"whipping" as the hull bends back and forth about the 
middle. With an explosion off the bow or stern, the shock 
wave hits the submarine as a plane wave due to the small 
cross-sectional area at the point of impact. Since the 
force of the shock wave iS axisymmetric about the centerline 
of the submarine, no whipping motion is induced in the hull. 
The inertial forces are transmitted down the hull until they 
reach the stern. The stern reacts violently to the combined 
inertial forces resulting from the explosion as well as the 
reflected pressure pulses from the water displaced by the 
translational motion of the submarine. After the initial 
transient response, the hull settles into an “accordion 
mode" about the middle where the motion of the bow opposes 
Piatto “the stern. | | 

The purpose of this investigation is to examine the 
shell/substructure interaction when subjected to UNDEX end- 
on loading using the ELSHOK code. The hull used in this 
study is a ring-stiffened cylindrical shell with endplates 
Similar in shape to those used as SSTV's. The internal 
equipment or substructure is a diaphragm attached to the 
shell whose thickness is increased in order to study how the 


mass ratio between substructure and shell affects the 


iD 


coupled transient response. The dynamic amplification of 
the response of the substructure at resonance 1S examined 


for both taper and conventional charges of equivalent 


impulse. 
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Il. EXPLOSIVE SHOCK RESPONSE USING ELSHOK 


A. GENERAL PRINCIPLES OF OPERATION 

The ELSHOK computer code consists of a family of 
programs developed to calculate the transient response of a 
submerged, ring-stiffened shell of revolution of finite 
length to an underwater shock wave emanating from an 
explosive source located at an arbitrary point away from the 
structure [Ref. 4]. The shell is assumed to be linearly 
elastic, with or without internal substructures, and the 
Surrounding fluid is treated as an infinite acoustic medium. 
Component modal analysis is employed in all phases of the 
calculations. The complete structural system is considered 
to consist of the ring-stiffened shell and any attached 
substructures. The-vibration modes of each component are 
calculated separately, and the equations of motion for the 
entire system are obtained by enforcing compatibility of 
deformation at the points of attachment. The free-free 
modes of the empty ring-stiffened shell and the fixed-base 
modes of each individual substructure are coupled through 
the use of dynamic boundary conditions [Refs. 5, 6]. This 
eliminates the need to calculate the combined modes and 
natural frequencies of the entire system as well as the 


requirement for a combined system stiffness matrix. 
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The structure-fluid interaction 1S approximated in 
ELSHOK using the Doubly Asymptotic Approximation (DAA) 
methods of Geers [Ref. 2], expressed in terms of functions 
which are orthogonal over the wet surface of the submerged 
shell. By matching exact sHescume= oleae relations at 
zero and infinite frequencies, the elements of the matrices 
in the resulting DAA are obtained. In transient problems, 
the DAA yields exact solutions at early and late times 
providing a smooth transition between these two limits. By 
accounting for the effects of the fluid with quantities 
defined solely on the wet surface of the shell, the DAA 
essentially uncouples the fluid field from the structural 


field. 


B. ORGANIZATION AND OPERATION OF ELSHOK 
As stated previously, the ELSHOK code contains a group 

of Computer eT Oanents which are utilized in order to obtain 
velocity “time histories for various locations in the shell- 
Substructure system. The programs, in order of execution, 
are: | 

1) BOSOR4--structural analyzer for shell [Ref. 7] 

2) ACESNID--virtual mass processor 

3)  PIFLASH--shell-fluid processor 

4) SAPIV--structural analyzer for substructure [Ref. 8] 

9)  PICRUST--substructure processor 


6) USLOB--time integration processor 
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7) PUSEOB=-pletting precessor 
Figure 1 illustrates the general relationship between the 
programs which make up the ELSHOK computer code. As can be 
seen, the operation is divided into four phases. 
1. Phase I--Shell and Fluid Analysis 

The first step of the ELSHOK calculation is to 
determine the in-vacuo free-free modes and natural frequen- 
cies of the shell. This is accomplished by modeling the 
shell using the BOSOR4 finite difference code. This code 
applies to segmented, ring-stiffened, branched shells of 
revolution having various meridional geometries, wall 
constructions, and ring reinforcements making it well suited 
for the modeling of submarines or SSTV's. For compatibility 
with ELSHOK, a separate BOSOR4 calculation must be performed 
for each Cir eum Sronten hetero i: (N) included in the 
analysis. Examples of different harmonics in the circum- 
ferential distributions include N = Oreo cdentnay cer sional 
and N = 1 (pure translation/whipping) as shown in Figure 2. 

The second part of Phase I is the computation of .the 
virtual mass array using ACESNID. This provides the 
late-time contribution of the DAA. The virtual mass array 
is determined from the solution, based on simple sources of 
a low-frequency steady-state problem in which normal 
displacements corresponding to surface expansion functions 


are applied to the surface of revolution in the infinite 
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fluid having the same shape and size as the net surface of 
the shell. ACESNID requires only one execution for all 
values of N considered. 

The execution of PIFLASH completes Phase I by 
combining the outputs from BOSOR4 and ACESNID to create a 
"shell-fluid file" which contains all the information 
required to describe the shell during subsequent operations. 
One of the useful features of ELSHOK is the fact that any 
number of substructures can be analyzed for a given shell 
without having to repeat Poaise lecal cit atiotor. 

2. Phase II--Substructure Analysis 
‘The analysis of the substructure is performed using 
the SAPIV finite element code. The SAPIV code contains a 
variety of finite elements (plates, beams, pipes, etc.) 
which can be used in the modeling process. It is a general 
purpose code based upon a lumped mass formulation applying 
to linearly elastic structures. The PICRUST code takes data 
from SAPIV and reorganizes it to facilitate the solution of 
the equations governing the transient response problem under 
Study. It 1s during the execution of PICRUST that the 
connectivity between the substructure and the shell 1s 
accounted for and influence coefficients for evaluating the 
forces developed at the points of attachment are calculated. 
3. Phase III--Submerged Shock Response 
In order to prepare a combined shell-substructure 


input for the USLOB time integration processor, the 


Ze 


shell-fluid file from Phase I is merged with the 
Substructure files from Phase II. USLOB is the portion of 
the ELSHOK code which contains the underwater shock from an 
arbitrary point. It allows pressSure-time history inputs for 
taper charge modeling or an exponentially decaying pressure 
impulse such as those encountered with conventional charges. 


henis Of the form [Ref. 9]: 


a) = x (wt/3/R) 2 exp(-t/@,) (1) 
where: 
P(t) = incident pressure on the shell .(psi) 
Ky = multiplicative constant for incident pressure 
Ko = spatial decay constant for incident pressure 
te = time after arrival of shock wave at point of 
interest (msec) 
W = weight of spherical charge (l1b.) 
Pw 2/3554 es , 
Oy = KW (w /R) = time constant of exponential 
decay (msec) 
K. = multiplicative constant for time constant 
Ky = spatial decay constant for time constant 
R = distance from the explosive to the point of 


interest (ft) 

Kit Ko, K3, and Ky are constants which depend on 
explosive type. USLOB employs a modified version of the 
Runge-Kutta integration method to produce velocity-—time 
histories in tabular form at user specified points on the 


shell and the substructure. 


ZS 


4, Phase IV--Plots of Velocity-Time Histories 
After the velocity-time histories have been 

obtained, PUSLOB is utilized to display the velocity 
responses in a plotted format on a TEKTRONIX graphics 
terminal. These plots are adequate for viewing trends and 
general responses, however for formal display purposes, the 
IBM Versatec plotter yields a superior product. Appendix A 
contains a code which was written to convert the velocity 
punch-card files created by PUSLOB into a data format which 
is acceptable for Versatic plotting using the EASYPLOT 


program. 
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III. MODEL USED IN THE ANALYSIS 


A. SHELL MODEL 

The shell used in this analysis is a small-scale model 
of a typical SSTV used for conducting underwater shock tests 
on high-impact Shipboard machinery. As can be seen in 
Figure 3, it is a high strength steel free-free ring 
stiffened cylinder with aluminum endplates. Appendix B 
includes a sample BOSOR4 input for one of the 
circumferential harmonic distributions. Data can be input 
to BOSOR4 using a formatted input code or by responding 
interactively to a set of prompts. Each segment is modeled 
separately which allows for different materials and physical 
properties to be included in the structure. The small 
stiffeners along the shell are represented by an orthotropic 
approximation which in effect increases the density of the 
shell. The six discrete rings are modeled individually, 
allowing for different size or material. The torsional 
rigidity (GJ) for each ring is found using the following 


relation: 


GJ = EJ/2(1 + v). (2) 


where: 


E = modulus of elasticity (psi) 


a 
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ei 
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at/(40 I) =—~soOrsilonale rigidmty constant: [Ref. 10] 

v = Poisson's ratio 

Since the shell is a surface of revolution, the 

properties of each ring cross-section need only be specified 
at one point. The end plates (segments one and three) are 
divided into ten equally spaced nodes, and the cylindrical 
portion of the shell (segment two) is divided into forty- 
five nodes along the longitudinal axis. As with any Anite 
difference code, greater accuracy will be obtained in the 
Somution with higher numbers of nodes, however, computa- 
tional time iAieineneasc. Siwo additional modes are 
automatically inserted by BOSOR4 into each segment in order 
to reduce the truncation errors associated with segment 
interfaces and to prevent spurious vibration Oe Due Cisring 
modes. Through finite difference techniques, BOSOR4 
determines the in-vacuo free-free modes and natural 
Paequencres OF the shell. Since this study deals strictly 
with end-on loading, the problem is axisymmetric about the 
longitudinal axis. Although only the N = 0 (breathing/ 
Seerolonaly modes are activated, a few N= 1 (translational) 
models must be calculated and retained due to the fact that 
the ELSHOK computer code was developed for side-on loading 
and, therefore, requires a translational input. Since there 
is no translation, the solution is not affected by these 


Superfluous modes. BOSOR4 generates a shell file which 
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gives the mode shapes for each natural frequency. All 
breathing modes have been retained for further study, 
however, the torsional modes have been dropped from the 
analysis since they are not activated by an end-on load. 
Table I contains a listing of modes and natural frequencies 
retained for input into the PIFLASH program. 

Figure 4 gives a representation of a typical shell- 
substructure configuration. The two coordinate systems used 
in an ELSHOK analysis are illustrated. xX, Y, and Z@ refer to 
the global coordinate system of the shell, while x, y, and z 
refer to the local coordinate system of the substructure. 

In both systems, the x-axes run longitudinally down the 
hull. The substructure z-axis is related to the global 
Z-axis by an angle «, where o refers to the substructure. 
Points on the wet surface of the shell are located by 
specifying the meridional arc length s, measured along the 
reference surface, and the circumferential coordinate angle 
8. The lower case u, v, and w are used for local shell 
displacements. To eliminate torsional modes from 
consideration, all N = 0 modes in the shell file with 
Significant v-displacements are discarded. 

Once the shell files have been obtained from BOSOR4, the 
input code for ACESNID can be created. This program 
determines the virtual mass array which produces the late- 


time contribution of the DAA. It considers a cavity in an 
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TABLE 1 


MODES AND NATURAL FREQUENCIES OF SHELL 
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FREQUENCY 


1. 259E=03 
2.2390E702 
3.850E+02 
4.7618+02 
Be Senor 
Feo LOE TO 
956 9H mo 
1.048E+03 


Ve OO 4 Eeu Ss 


1.585E+03 


eo? Ont 


1.598E+03— 


We OG Ors 
ie sae OS 
1. Oke + Os 
ee Orerornere Ome 


1.878E+03 


ore 0 eee 


io Gees 


1.390E-04 
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imeinute acoustic fluid corresponding to the wet surface of 
the shell. A user specified number of surface expansion 
HuneEtons heying the property of orthogonality over the wet 
surface of the shell are obtained during the calculation of 
the virtual or entrained mass. The cavity is divided into 
bands which describe the behavior of the fluid. A 
sufficient number of bands must be provided to enable the 
fluid response to match — normal motion specified by the 
surface expansion functions. In a sense, they can be looked 
upon as nodes in a finite element model. In this analysis, 
the shell is assumed to be totally immersed in salt water 
with a mass density of 9.59684E-05 lbf-s*/in’. The 
endplates are divided into fifteen bands each, and the 
panacea: has fifty-one bands. Six surface expansion 
functions have been generated for the endplates with ten 
along the cylinder. | 

The shell mode files from BOSOR4 and the virtual mass 
file from ACESNID are combined in the PIFLASH program to 
produce a shell-fluid file for further processing. The 
model of the shell is now complete and can be used for any 


variety of substructures without any recalculation. 


Dee oUbo RUGTURE MODEL 
The idealized internal equipment selected for study in 
this analysis is a high strength steel diaphragm or plate 


located at the discrete ring at frame nine. The SAPIV 


al 


finite element code is utilized to model the diaphragm. The 
diaphragm has been divided into eight "pie-shaped" wedges 
which are further subdivided into three sections. Twenty- 
four plate elements are used in this model with twenty-five 
nodal points. Since this study is concerned with end-on 
loading, the node at the center of the plate is constrained 
to move only in the x-direction with no rotations. The 
outer nodes Wee through twenty-five) are rigidly 
attached to the shell. This is accomplished through a 
modification of the general purpose SAPIV code by the 
developers of ELSHOK. During the parametric study, the 
diaphragm thickness is uniformly varied from .25 inches to 
9.25 inches. Figure 5 is a TEKTRONIX representation of the 
finite element model of the substructure. Appendix B 
includes one of the SAPIV input codes used during the 


parametric study. 
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IV. ANALYSIS 


The purpose of this analysis is to study the shell/ 
substructure interaction after the system is subjected to an 
end-on load while varying the mass and stiffness of the 
internal substructure. Two different types of loading are 
explored: a taper charge and a conventional charge. Taper 
charges are characterized by sustained incident pressure on 
a body over a period of time. They are used in underwater 
shock tests to simulate the type of pressure profile which a 
body would experience from a nuclear detonation. Extensive 
use of scaling is involved in efforts ae opredice the 
physical’ behavior of an Obj]\ece tamrcspousec =o anc lear 
charge without generating the extremely high pressures 
associated with actual detonations. Conventional charges 
create a large initial pressure pulse followed by a rapidly 
decaying exponential. ELSHOK uses inputs of charge weight 
and standoff geometry to calculate the transient velocities 


of the shell and substructure in response to these loads. 


A. ANALYSIS PROCEDURE 

Once the shell model has been produced, a series of 
SAPIV runs are conducted for the diaphragms of increasing 
thickness. Each run generates natural frequencies and 


accompanying mode shapes for each model. The PICRUST 


34 


program reduces the data obtained through SAPIV into a 
format which can be combined with the shell-fluid file for 
processing with the USLOB code. The Bipstederure 1s 
connected to the shell during this phase of the ELSHOK 
analysis. The diaphragm is attached to shell segment two at 
frame nine (node eleven). 

The weight, type, and location of the charge is 
specified in the USLOB input code. To set up an end-on load 
referenced to the global coordinate system, the charges are 
placed in the negative X-direction at a distance of 70 feet 
(840 inches). This distance is typical of stand-off 
distances used during actual tests. For the taper charge, a 
typical incident pressure-time history shown in Figure 6 is 
used to simulate the shock wave loading. In order to make a 
meaningful comparison between the structural responses from 
a taper ataree and a conventional charge, it is necessary 
that they have equivalent shock wave impulses. In the field 
of underwater shock analysis, "impulse" refers to the time- 
integral of the pressure profile [Ref. 9]. The impulse of 
unit area of the shock wave front up to a time t after its 


aaa alelbs Givens by: 
ie 


I(t) = | 2 (a) felis (3) 
0 


The impulse of the taper charge is found by calculating the 


area under the pressure-time history in Figure 6. AS was 
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Figure 6. Incident Pressure-Time History--Taper Charge 
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stated earlier, a conventional charge follows a pressure- 
time history given by equation (1): 


ao 


1/3 7R) * exp(-t/®,) (1) 


P(t) = K, (w 


where: 


2. Kawl/3(wl/3/p) 4 = time constant of exponential 
decay (msec) 

HBX-l1 1S a popular explosive used in underwater shock tests, 

SO it was decided to use it in this study for the 

conventional charge. Ky Kor Ky and Ky are constants which 

depend on the type of explosive used. The values used in 


thas study wause wbRef. del]: 


K, = 3.8354314E+05 
K, = 1.144 

K, = 3.03131E-05 
K, = 7.247 


By matching the impulses between the taper charge and the 
conventional charge, a charge weight of 352 lbs of HBX-1 is 
found to yield an equivalent impulse. The peak pressure 
generated by the charge is 1620 psi which creates the 
incident pressure-time history shown in Figure 7. The 


conventional charge expends most of its impulse in the first 


ou 


millisecond, while the taper charge takes over four milli- 
seconds to expend an equivalent amount of impulse. Because 
of the sharp rise in pressure experienced from a conven- 
tional charge, a higher transient response is developed in 
the shell than with a taper charge of equivalent impulse. 

After selecting suitable charges, the time step incre- 
ment and integration limits are specified in the USLOB input 
code. Through a trial and error process, enough time steps 
are chosen to identify the significant interactions between 
the shell and substructure. It was found that 1600 time 
steps for an 80 msec time period gave an adequate illustra- 
tion of the system response. The USLOB program enables the 
Neen des ot anene the velocity-time history Meeponse of 
various nodes in tabular form. This gives some indication 
of whether the run was successful, however, the plots from 
PUSLOB are required before response characteristics can be 
identified. 

To obtain a baseline velocity-time history, a complete 
analysis is performed on the submerged empty stiffened 
shell. From the velocity-time history plots, the dominant 
frequency from the computed shock response is observed. 
Provided the mass of the diaphragm remained small with 
respect to the shell, it was felt that this provided a good 
estimation of the excitation frequency for the diaphragm. 


The dominant frequency is 187.80 Hz for the taper charge and 


3e3 


195.44 Hz for the conventional charge. A series of SAPIV 
calculations are conducted on diaphragms of varying 
thickness in order to identify the fundamental frequency of 
each substructure. These values are given in Table II. An 
analysis is performed on the shell/substructure system for 
each diaphragm thickness to observe the coupled velocity- 
time history response as the substructure mass grows in 


relation to that of the shell. 


TABLE II 


FUNDAMENTAL FREQUENCIES OF DIAPHRAGM SUBSTRUCTURE 


Thickness (In) 
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Frequency (HZ) 
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V. RESULTS 


The velocity time-history response begins at time zero 
when the incident pressure wave from the underwater 
explosion strikes the shell. With end-on loading, the 
forward endplate receives the initial contact with velocity 
indicated along the longitudinal axis. The empty shell 
response is analyzed first, followed by a discussion of the 


coupled effects when the diaphragm achieves resonance. 


A. VELOCITY-TIME HISTORY RESPONSE--EMPTY SHELL 

The empty shell velocity-time history responses for the 
conventional and taper charges are contained in Figures 8 
through 20. It can be seen that the conventional charge 
causes a much greater transient response in the shell. In 
both cases, the velocity is rapidly damped by interaction 
between the shell and the infinite fluid medium. It is 
interesting to observe the differences in response between 
the various locations on the shell. The forward snaomeee is 
perpendicular to the shock wave, receiving the full force 
from the incident pressure pulse. Inertial forces are 
generated by the endplate causing severe local deformation 
in the nearest nodes on the cylinder due to the rigid 
coupling between segments. The motion propagates down the 


shell along the longitudinal axis until it reaches the aft 


41 


zezUeD ezeTdpug qMd--(TTeus A3duyq) esuodssey AzTORSTH SUTZ-ARTOOTAaA *g oanbty 


(ONS) SUINEL 
og OZ 0g OS OP OE . O02 oT 


S 





Or- 0S- 


Oc- 


Rules 8)¢. 8 6), see ENO 6.6 G9 08 6\o.e1eim 6 (Oe Meleren *)ee radia (6 50° SUIS lee) S00) g enere 
ME CRIVTGTeTG:DIUls) LeuelOlekelelo Ol S16.*) WO lisiei~'9/619.6:-\018 8:4 lwiessie:« s\eiateletem ee = 10 /0'S 0.610 ni @ieimme Ue  aguesw se) ene 0 Ce Sere Soe ° 8.816 © Rie ee ime Sees So ee 


0c- 


TA 





U 
(ORS/LA) ALIOOT 


OT- 


of 


0¢ 





. 
Cee OE SE A ew eee ee sseeeeeerteeee 


- 


0c 


» 
=”. 
eoce 


GOUVHO Wad VL 
GQNd9gT 





ANO GON — ANO LNAWDAS 
TIAHS ALdWa — ASNOdSHa ALDOTHA 


42 


OlezZ owerig--(TTeus Aqjduy) sesuodssey AAOASTH sWTL_-ARTOOTSA *6 aanbty 


(OASW) FACL 
On a2 Og OS Ov Ge O02 






seer eens corerreee «= 


OP eer cece n tance aemeogtoneorrteeerereseererrs-Oere He en 8 


cI~ 


Oi 


(94S/L4) ALIQOTA 


¢- 


i 
v 


G 


a 
MDVD TVNOLLNGANOD 8 
qOYVEO BEdVL 
ANSE 
n 


OMS ACON ~ OAL LNANDSS 
TTHHS ALUN — ASNOdSHY ALTOONA 


43 


SATY owerq—--(tTteus Aqduq) esuodssey AAOISTH BUT L-AFTAOTSA 


(OASW) FALL 
og OZ og OS OF o£ Uc Gt 


_— ASUVHO TYNOLLNAANOOD _ 
GOUVHO Ydd VL 
QNYONT 


NAHAS GQON — OAL LNAWDdS 


TIHHS ALGWH — ASNOdSAY ALISOTHA 





in 00h oat 


o> OT— 


0 
(QHS/LA) ALIDOIA 


44 


SUTN SUeATT--(TTeuS Ajdug) oesuodssy AAOASTH suTZ-AqTOOTSeA “TT sanbty 


(OASW) ANIL 


08 OL 09 09g OF O£ 02 Ot 0 
Ma | hf 
Z Bc “ : 03 
: ; : on 
ane Pnee coePnen stig atalevetars sta aleis ee aie ereimmein elec eieves pstersls ee Stee ee nei ie wlelsl clateleieereislcleleic/ale/eeleielevetelel iglcinmtsitieicisis/ev/®!e.= 01) Mi n@ueisiteceiniacc) ie) Uc islelclelel #o\eie(ikeMmel@meleicioue eile nO 
Bee atte het e ; : ny 
: : : | 
Dp ef stata ects Bale ahs Rte Tics Ae space Haan a eee ee oe aia ashe nem etits # RIOT te eee ales ecient eet ale ro ere eee reliome eae ae 
ene reheat eulaaas 7 
: 
SOEPRE He) Had ca deh vies crensneewes) «ulerteusenens siseesile SusigeGina<s ic useu@setedescesrwssewedincssgas weeant+t0es ra) + sgeqmernen ie \ seas eaauiaan ccs = 
: n «<4 
: a 
; nN 
‘ ve c= 
d Sic Ph MRS. 2 «tise. ane Rea a co Ree Ee aa eee is antec Fe eee i© 
oo ae te eee Seb ae Cuaee a eaten weneses cus. “Saauee - an = 
' mi 
' 
a 
' ~< 
: on, 
' 4 
‘ 
' > 
UV) 
t. 
' tz4 
1 ©) 
' Seen” 
4 
4 


QNA9D4T 





NAAUTA AGON — OAL LNAWDS 
TTHHS ALAN — ASNOdSHY ALIOOTHA 


45 


UsesRITYL swearajg--(TTeus Ajdugq) esuodssy AAOASTH suUTL-ARTOOTIAA 


(OaSW) ANLL 
og OZ 09 OS OF O€ O02 OT 


—_ ADUYHO TVYNOLLNGANOOD _ 
AOUVHO Ydd VE 
GNGDET 


NaaLald AGON — OAL LNANDES 
TIAHS ALKHNA — ASNOdSHa ALDOTHA 





"ZT aanbty 


Oc- Sce- 


i 


Ol= Sie 


S- 


AS/LA) ALIOOTHIA 


o 


( 


46 


usezybtg owerq--({TTeus Aqjdugq) esuodssey AzoAsTH auTZ-ARTOOTAeA “ET eanbhty 


(OMSK) FALL 


og 02 a9 09g OF O€ 02 ol 0 
: : i] 
: : : : an 
BT aie oven oer e W Alaieis Sistine elicieve so 0 sc(@jerp 0! Blais 8 Sieteaie alcies @ereieicrs0l-lor hs ieelieer= esis Se eee eee tales ae loleineee Cee OO cere bee wie ce cte acy oease eeV@im lore eo eneleraie ee aoe exe n 
: : : ‘ © 
3 : ee eee ae ! 
Bee pager OO Ia Saale” le we OOO Oe Bee Gee ee Cer eee ee rob rats earns. - 
3 : ! 
: ~ 
Gu 
: tJ 
a SS 
BT avale ete cruciate late lolnelbletecsctets lave eiaissacoiciare es8 Oia Bis alerslo eco 1861816 a.6 818 jctthalsicvers e8@e e@8e..+4e 0 8e 8+ oO 8. es he tes nue @ e-. CGete ws oe Bs ees aan as a oP ae ae See Gl ve wee gn OQ 
pean 
J 
me 
— 
4 
Bu 
W* 
trJ 
a) 
“Nee” 


_ AOYVHO TYNOLENTANOD _ 
GOUVHO Wad VL 
aQNqod 





ALNAM.L GON — OAL LNAWDIS 
TTHHS ALANA — ASNOdSde ALDOTHA 


477 


seayp,—-AquemMy, ewerzg--({TTeus Aqduy) esuodssy AAOASTH suUTZ_-AQTOOTSA “PT sanbty 


(OUSW) AWLL 
08 oz 09 OS OF OF 4 OT 0 





a ema oeeemeeremeacanae es t 
“ : se 5 9 
: : : : ; on 
: : : 
ee emer e rete > cc ele mn otncvoleiotetcleleierou lm viliselcfatel=t-Vel=(elels) =/7l-/=1aiiislelujs\w\slexallitelis/e/u/eliei®/e) | (tel -/-1e/-N-/@\-Le1-1°1°)9) --bel-feiesesniainLels|¢ ain oe s/~ieiateelcl= Sl rn 
° ° . Oo 
{ 
Oe Qeer enon eeer "ess on ee er rr ce i ete mt tee Oe See or ener ots  OFHH SEF OOE . erreeeuce POSE SOOO OOS erent eernes bom 
. nn 
: J 
aI cs SIRE IR TE Te reiei eis reai mea areca ee enter eve Meanie OMS oterne tclerecslnieloaloreinielerstiemmsters|ercicinis cieiet (<M SteTe (eer tacetcas “ Pee ee 
: So <i 
: ie 
| 1o 
2d DOS DERE BCH CL 0016 CUDA CORD EOOnDe Ses HACOCDaSCOCne mare : on S 
a oat oa 
ts | yj 
’ 
a 13 : 1 oO ~< 
Preraveicicve) <a virat ~ciaietelntotalelmletsovencierstele elerecselsreais/crermis/alsieis) sisi nniel Wucicleleteterei@iel sic (soi6lels/eisiie\=(eis//aVal lw alei==10)0 =0)n (01597 0Il s\lp)ro10 1616 nie leiei sure oe ; ve 
‘ : ‘ x] 
‘ 
e ' i ‘ — 
' = 
: : ! Y} 
Mipi@eiminlere (ats) in) =)= =) =) =in\=/pi8160n)6)e[al |= lao /e 601 ® = @inis Precise Sse es a plaieleiere aie @ie/t uleier cue l® lel eit@teleiahsc wiv) wla/s’m) (se falere "Gly mio 016m @.0 (0) O/luia of) oe! pees ie z meteors : nn by 
: ane ' © 
s : A r? 
5 ‘y ,! Pee 
tipieisteleisiii=is(0\010/4)s/slesielole sle)siciele|o/ala|sie's ele'e six sWreieiavelctateleferelorelsleieeitctcie(oisle Mele =(e sl elelslelel=1-Velsieteit-ictal=la1e]=i telson ielel oisi-/siels)olel=1-1elo(s/a 16/6 (=a) @iel=l tals /sle = iinisis/s(ei=1e\8(e) noon ni ‘; eneg Oo 
: y i! 
: a 
if 
— 
“ON 





GHOUVHO Wad VL 
GQNdDaT 





GLLNAANOD | , ke ee 





YNOA ALNAML AGON — OAL LNAWDaS 
TTHHS ALdWa -— ASNOdSAY ALIOTHA 


48 


usAdss-AquomMy, swerj—--(TTeus Aqjduy) ssuodsay AAOYSTH swtpZ-AAtoojTsapn 


(O4SW) ANIL 


og OL 09 0G OF OF 02 Oy 0 


GAOUVEO ddd VL 
(NAOT 


LHI ALNIML — OAL LNIW9AS 
TTHHS ALdWa — ASNOdSHaY ALDOTHA 





of S 0 cS-at- 
(O2S/LI) ALIOOTHA 


ST 


GT eanbty 


49 


OM,-ARATYUL eWeAg--([TTeusS Aqjdugq) esuodssy AAOASTH swTL-ARTOOTOA *9T SANbtTy 


(OUSW) ANIL 





0g OL 09 OS OF O£ 02 OF 0 
4 _ 
: an 
! 
aM ood ocr oc, nee eared acid cae. o8:0ike's dn dR oo eRe a CTS a a c 
oO 
Be scsi delidiciic: Mele CEE sone ae.n: RM MEMUENE  tosidhurdwas cucu 4 Um ele saat a ee ul 
WE en een occa ences SO ag ee ea siaice ee 7 
re: eee. a. ne iu 
I eee ee oy ee aoe GE ep eanG) alnn do SCCM eR Se swe pa oe EMM ain oh oA SRG oe ere 5 
: ; tr] 
: - a 
I ee ee eee’ g hehe ees shece sets as ogee Bacto lye chs iso a elect as ca olga 5 clas RE MRR cs nn aN Age 
ae ve 
‘ fc ‘ ~j 
ya 2 
Ve ‘ or. 
| = 
: 3 
_ 0 8 
BD 
° 
O 
° 
AVUVHO Wadd VL 
QNGSaT 
nN 
on 


AAMHL ALULLL 3GON — OAL LNENDS 
TTHHS ALANA — ASNOdSaHe ALDOTHA 


50 


XTS-AQATYUL SWeIgq—--(TTeus Aqduwyq) osuodsay AAOASTH SWTL-ARZTOOTOA 


(OASW) TNLL 


0g 02 09 09 OV OL 02 ot 0 
‘ 
a nN 
: am 
| 
Marre) icy ctavalavelotaievelsvalelelsiavelarslajsielets alsierer) 6 oie)» el eferave\mcsie Wal scieceneleie lars) 6: © iay®\n mevelSisimn-K@om | oaucisyonere\® ie) 816) $) (6 8i~(@/ei/e\uiaiiw\cia = 5%] °/\r! sissy OS 512.900 aaa aca ac n 
© 
ft 
We aca a ofa a diss nig dvaia is or0:6 6 0.8 aera eS Te ae tiv) Ee ene aie seiela rie welcieils SCAG iginibig-rfeiairis me wen ciereicine cia?! > nghesie cia ne miler mis eciei8 cleat eres) ste eisai ie 
om 
1 
ee ere io iniictey oy aereiele eratetere re iavels ave elated olers ole clePmelsPacsisia eve volar late oimctitelstetale/svete eie1orsleirMoulvrelcyeroiKe\siorsleje siete eieYoiniey©ieiskejslerenmateL Tees uese 7 = : ae 
° o 
.: 
e i 
au i 
is oo 1 
alalataleis aie sisi lclelel:| wie eielsrele\e) sis 00) 6 oe eelele lea wing) aise 6s\@ 6.61016 .a.ujs) O16 2's) = ola 10 0la ee Olsens (maui cex0le,= a) 010 0.18 @) )10)8),0/616 010 6:66) 0/ Oi) ie 8s8.01e:8.e/6106is/ ee ieee ae t a oh on 
‘ ‘ 4 
‘ i: to 
t t $a 
‘ a | 
4 x t 
‘ t 
‘ t 
’ 
’ 
t 


_ADUVHO TYNOLLNAANO) _ 
GOUVHO Yad VL 
GQNaoa'T 





NAATS ALUIHL AGON — OAL LNAWDAS 
TIAHS ALAGNG — ASNOdSHa ALDOTHA 


(OaS/LA) ALIDOTHA 


"£T vanbty 


ay 


Azz0qg owerg—--({TTeus Aqdwy) esuodsey Az704STH SUT L-ATOOTSA "QT aanbty 


(OAS) TALL 








0g OZ 09 0S OF O€ 02 of 0 
: z “1 
on 
{ 
RN AM oc eM OC oe “a SUM et we'd o si Gee ae Sk 0 9a toe Se ae aeRO, oA 69 
°o 
ssain'we’necud tS RO Lao Sa PR al 
TR 2 a |, oS PME vs ing sloglti tama Bree eR ee Oa eT Rm he die wate ca o 
1 ee eee to 
' on 
\ rr 
5 
es a 5 
\ wl 
! < 
i o 
ee 
: 4 ‘ ~~ 
: inf by 
: a } trJ 
: aia ‘2 
: ioe ‘ — 
; | 
: rf | = 
: ‘ 
: i 
: a 
: t baa 
: Mi 
" 
‘, 
; no 
oe 
GND 
nD 
On 


ANO ALUYOd JOON ~ OAL LNAWDTS 
TIAHS ALdMNA —- ASNOdSHY ALIOTHA 


Sy 


aatq-Aqaoq owerq--(T Teus Aqdwyq) esuodsey AAOISTH QUT L-ARTOOTSA 


08 


(O4S) GINIL 


04 09 as Or O& 02g oT 0 


_ SOYUVHD STIYNOLLNAANOD _ 
GADUVHO Wad VL 
GNA aT 





KIS ALUOd AGON — OAL LNAWDaS 
TTWHS ALdWa — ASNOdSAY ALDOTHA 


oz st oF S 0 COT 
(O4S/LY) ALIQOTHA 


Se 


"6T oanbty 


53 


qequeD ereTdpug Lav--(tteus Aqdug) osuodsey AAOASTH SWTL-AQRTOOTSA *0Z oanbty 


(OASW) AWLL ; 








og O24 09 0S OF Of Od ol 0 
: uw o 
© a 
: vt 
: as . 
reterntta icles’ | citinieiy) 10 6is:wlelareinininiersls isle 0 /e'e'> 018 @i-1-(e1eele\si 16/010 61n)8 n\@n\© .Oin'ea!0is\a 6 nioinls 01010, @\0 @:8)oi8'S 8.101 41+ 316.610) -.0(8)¥'9/0.8)s ig Als) a -eleinielsieisicieys) «isle isis Oc elele ei eis io ied’ oie eis) cere ee Sie 00 sOrerece a 
: 1 2° 
2 4% 
: es 
= 4 j 
Rohe mere HPCE OOHRS HEN OH SESE SEH ete Hr FFP OHL ees SOEEO® POC ere eee ere reeset FB |r Se OH Her orotate He gs "Ft erersserrere® sheicieieieisie sis /eitieeisieivisisis eooore ¢4 é3 
; oe oS 
: at 
‘ a4 
: i 
Pee: ney Oe - 
Mae aan O 
ie : 4 4 te 
“oes aa rites 
Vice. MN hegre pon © 
‘: i on 
1 re yt — 
1 bee 4 3 
' : ThA ~< 
t t a -9O L~ 
ieevet ae euler 7 
eo Oea <a 
tee ae - WY 
i : « fe : o 
$ >. oUt OW de 2) 
¢ a Wane Ty — 
vt hd! 
; ‘ pale eistaiclere 2 ; = 
4 \y ' 
| 1, 
u _ 
y m 
¥ 
ee ee Iss su dees coker ee 
Aq)UVHO yaaVL : tf 
GNA : if 
: : : 4 on. 
Qo 


AATAAL AGON — F4UHL ENEWDAS 
TIHHS ALdGWY — ASNOdSHaU ALDIDOTHA 


54 


endplate over a millisecond after initial shock wave 
contact. At first the aft endplate moves in the negative 
X-direction as the "ripple" proceeds down the structure. 
Then it experiences a violent transient response as it is 
affected not only by the inertial forces from the shell, but 
by the reflected pressure pulses from the water it 
displaces. After the initial shock, the shell settles into 
an "accordion" effect, expanding and contracting as the 
forward and aft endplates move in opposite directions. 
B. VELOCITY-TIME HISTORY RESPONSE--COUPLED SHELL/ 

SUBSTRUCTURE SYSTEM 

From the results of the analyses performed on the coupled 
shell/substructure system, Tables III and IV were created in 
order to identify when the diaphragm achieved resonance. 2/2 
is the frequency ratio between the fundamental frequency of 
the simply-supported diaphragm (2) and the dominant frequency 
observed in the velocity-time history response of the sub- 
merged empty stiffened shell (2). M1/M2 is the mass ratio 
between the substructure (M,) and the shell with endplates 
ieee toune 21 1s a plot of the velocity of the diaphragm 
center versus the frequency ratio 2/2. It reveals that in 
both the taper and conventional charge responses, the 
diaphragm achieves its maximum velocity or resonance when the 
fundamental frequency of the substructure is nearly equal to 
the excitation frequency from the shell. Resonance occurs 


when the mass ratio is approximately .08. Since 
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MAXIMUM VELOCITY RESPONSE OF DIAPHRAGM SUBSTRUCTURE 


Thickness 


o25 


90 


TABLE 


len 


WITH VARYING THICKNESS (TAPER CHARGE) 


(In) 


2/2 


2a 
.465 
Gre 
930 
.976 
986 
995 
loge 
1.069 
1.162 
2.092 
3 e022 
3.952 
4.881 
5.811 
6.741 
7.670 


Seal 
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M1/M2 


aopite 
036 
0515 
.073 
.076 
S077 
.078 
080 
.084 
091 
.164 
1237 
s Sako 
383 
456 
.528 
.601 


~674 


Velocity (ft/sec) 


22.1.9 
16.64 
20S 
34.81 
3 SOEs 
Sec 
3 o> 
33S 
3.0i.310 
26.24 
1 dail: 
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TABLE IV 


MAXIMUM VELOCITY RESPONSE OF DIAPHRAGM SUBSTRUCTURE 
WITH VARYING THICKNESS (CONVENTIONAL CHARGE) 
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Thickness (In) 2/2 M1/M2 Velocity (ft/sec) 
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STIFFNESS RATIO (K4/K2) 


Stiffness ‘Ratio for Various 


2-DOF System First Modal Frequency vs. 
ih 


Mass Ratios [Ref. 


Figure 23. 


f= @ =A KI/ME = 3/7 RO/M2 (4) 


where: Kl overall stiffness of the substructure 


K2 


overall stiffness of the shell 


1t can be said that: 


K1/M1 K2/M2 = K1i/.08 M2 


Or 


K1/K2 08 


The coupled shell and substructure can then be modeled as a 
two degrees of freedom system as shown in Figure 22 where Ml 
and Kl represent the substructure and M2 and K2 represent 
the shell. The large mass (M2) has an uncoupled natural . 
frequency (f.) which is equal to rf K2/M2- This is 
equivalent to ereyceniname frequency (%) of the empty shell. 
Since this 1S a two degrees of freedom system, it will 
possess two natural frequencies for the coupled system. {[t 
has been found in coupled systems that for decreasing mass 
ratio (M1/M2) and for increasing stiffness ratio (K1/K2), 
the first and second modal frequencies of a two degrees of 
freedom system approach the decoupled frequencies f and f. 
of two single degrees of freedom systems, respectively [Ref. 
12]. The first natural frequency is bounded by for and fy 
1s bounded by the second natural frequency. From Figure 23, 


it is seen that with M1/M2 = K1/K2 = .08, the dominant 


6l 


frequency of the shell 1s approximately equal to the first 
modal frequency (f.4) of the system. The result is that the 
shell/substructure system achieves resonance when excited at 
the dominant frequency of the empty shell. 

With the diaphragm at resonance, Figures 24 and 25 
illustrate that the transient response of the diaphragm is 
much more pronounced than that of the shell due to the large 
differences between the two masses as well as the soft 
spring between them. Initially the diaphragm 1s drawn in 
the negative Re chieecaion by inertial forces as the shell is 
displaced by the shock wave. With the increased motion of 
the shell at frame nine, the diaphragm receives more and 
more energy. At approximately t = 20 msec, the diaphragm is 
at its maximum velocity while the shell is relatively calm. 
At this point, the diaphragm begins to transfer energy back 
to the shell. Since neta energy is equal to 1/2 MV, the 
diaphragm is a small mass with high velocity while the much 
larger mass of the shell moves at a lower velocity. The 
energy of the system is transferred between the shell and 
Substructure as the overall motion is damped by interaction 
between the shell and the fluid medium. The transfer of 
energy between the shell and the diaphragm creates an effect 
known as the "beating" phenomenon. This effect is seen in 
the velocity time history responses for the case of 


diaphragm resonance in Figures 26 through 39. The effect is 
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quite pronounced at the endplates, because they act as "hard 
spots" receiving the energy transmitted through the shell. 
The cylindrical portion of the shell acts as a conduit for 
the energy transfer. Again, the nodes on the shell nearest 
the endplate experience some local deformation due to the 
rigid coupling between segments. 

The empty shell velocity-time history response 1S 
compared with the resonant diaphragm case in Figures 40 
tmeough 45. Of particular interest is the fact that the 
Vemocltyapmsofiles are virtually identical for the first 
fifteen milliseconds. During this period of time, the shell 
1s receiving energy from the shock wave and transferring a 
portion of it to the diaphragm. Because of the differences 
in mass between the shell and substructure, not much energy 
is lost from the shell. The peak velocities do not vary by 
more than five feet per second during this time period. The 
difference arises at t = 35 msec when the shell receives 
energy from the resonating plete. The decaying velocity 
response is disturbed during the exchange of energy which is 
observed in the beating effect. 

As the thickness of the diaphragm is increased beyond 
the resonant condition, the interactions between the shell 
and substructure are reduced dramatically. The shell/ 
Substructure velocity-time history remains within 10% of the 


empty shell response until the diaphragm thickness is 2.25 
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inches thick or M1/M2 = .164. By the time it is 9.25 inches 
thick, it represents over half the mass of the shell, and it 
acts as if it was an integral part of the shell slaved to 


the motions of the forward endplate. 


8 6 


VI. CONCLUSIONS 


Upon inspection of the velocity-time history responses 
at various locations along the shell, one can see why the 
potential for damage is so great from an underwater 
explosion off the bow or stern of a submarine. While the 
cylindrical portion for the most part serves to transmit 
inertial forces with little velocity response, the 
extremities of the hull experience violent transient 
responses which may be passed to attached internal equipment 
located in those areas. The inertial influence of the mass 
Ont — endplates on neighboring nodes of the shell is 
apparent in both the empty shell and resonant cases. The 
response from a conventional charge is seen to be much more 
peenouncéd han cthat of a taper charge of equivalent impulse 
due to ene rapid expenditure of energy against the forward 
endplate. 

The results obtained from the resonant diaphragm case 
Support the theory that with decreasing mass ratio between 
Substructure and shell, the first and second modal 
frequencies of the two degrees of freedom system approach 
tne decoupled frequencies of the substructure and shell in 
two single degree of freedom systems. This allows the 


motion of hull to drive the diaphragm into resonance when 


yf 


the dominant frequency of the shell motion is roughly equal 
to the fundamental frequency of the substructure. It is 
possible that this could occur during an underwater shock 
test of a small piece of equipment ina large SSTV. This 
might result in the failure of the equipment in a scenario 
which would not occur in a real submarine with its myriad of 
internal masses and substructures. 

The large difference in mass between the substructure 
and the shell results in a tremendous dynamic amplification 
in the response of the diaphragm at resonance. As energy is 
transferred between the two bodies, the large mass of the 
Shell causes its velocity to be quite small compared to that 
of the diaphragm. As is seen in Figures 40 through 45, the 
initial resonant case velocity response is virtually 
unchanged from the empty shell response. If the mass of the 
internal equipment was sufficiently small in comparison to 
an SSTV used in a series of shock tests, one ELSHOK calcu- 
lation could be performed on the empty shell to generate a 
set of velocity gospancé curves. The velocity-time history 
response at the attachment point could then be input into 
the SAPIV code for the finite element model of any number of 
small pieces of equipment, and the substructure transient 
response could be obtained without having to run a separate 
ELSHOK calculation for each shell/substructure system. This 


would result in considerable savings in computer time. 
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Although ELSHOK is a powerful tool which can be used to 
accurately predict the results of underwater explosion tests 
on equipment, some shortfalls exist in the program which 
limit its use in design applications. At the present time, 
there 1S no direct method to take the effects of the surface 
cut-off wave or bottom reflection into account. The surface 
cut-off wave forms in shallow underwater explosions when the 
compressive shock wave reflects off of the surface of the 
water cauSing a tensile wave to propagate. Bottom 
reflections form in deep explosions as the shock wave is 
reflected off of the floor of the ocean. ELSHOK assumes 
that the explosion is occurring in an infinite acoustic 
medium with a Single shock wave. Another weakness in the 
program is the fact ee cannot process oblique 
explosions; charges must be placed in a pure end-on or 
side-on configuration. Shock tests conducted in accordance 
with MIL-S-901D are either end-on or side-on loading, 
however, so this is not a serious drawback in predicting the 
transient motion of equipment in SSTV's. 

It is recommended that further study be conducted in 
this area using more refined shell models resembling full- 
scale submarines. Once a shell model is produced for a 
given submarine class, researchers can create different 
finite element models of actual equipment to observe the 


interactions between the submarine hull and equipment of 
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varying mass and shape when subjected to end-on or side-on 
loads. It is felt that ELSHOK is an effective method to 
give an indication of the response of equipment in 
underwater explosion tests without needless destruction from 
unexpected failures. Its continued use prior to UNDEX tests 


will lead to savings in both time and effort. 
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APPENDIX A 
PROGRAM USED TO CONVERT PUSLOB PUNCH CARD FILES 
TO FORMAT ACCEPTANCE BY EASYPLOT 

Along with the TEKTRONIX plots of the velocity-time 
histories, PUSLOB generates velocity punch card files for 
each requested node. The punch card file consists of a 
block of velocities in exponential format the size of which 
depends upon the number of time points requested. The 
EASYPLOT program developed for the IBM 3033 by Mr. John 
Mainwaring is a quick interactive method to produce graphs 
using DISSPLA without having to create a plotting program. 
It will accept up to 100 data points in tabular form 
provided it is in a column format for X-Y input. The 
PLOTCONV1L code was created to accept data from two separate 
pumeh card files Pe taper and conventional charge 
comparisons and combine them into one data file for each 
node for input into EASYPLOT. 

In order to prevent data overflow, the time step must be 
sized in USLOB and PUSLOB to generate no more than 100 
velocity data points. Both the PLOTCONV1 VAX 11/780 VMS 
computer command file and the FORTRAN code are included in 
this appendix. The command file must be modified prior to 
each run to specify the punch card files being read and the 


titles of the new data files. 
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NANAAN 


PLOTCONV1.COM - - COMMAND FILE FOR EXECUTION OF PLOTCONV1 


> 

3! 

9 ASSIGN PUSL4&.PUN FOROLL 

$ ASSIGN PUSL49.PUN FORO12 

$ ASSIGN MAXSHL1.DAT FORO21 

$ ASSIGN MAXSHL2.DAT FORO22 
9 ASSIGN MAXSHL3.DAT FORO23 
9 ASSIGN MAXSHL4.DAT FOQRO024 

9 ASSIGN MAXSHLS5.DAT FOQRO25 
9 ASSIGN MAXSHL6.DAT FORO26 

9 ASSIGN MAXSHL7.DAT FOR027 

9 ASSIGN MAXSHL8.DAT FORO028 

$ ASSIGN MAXSHL9.DAT FORO29 

$ ASSIGN MAXSHL10.DAT FORO30 
$ ASSIGN MAXSHL11.DAT FORO31 
$ ASSIGN MAXSHL12.DAT FORO32 
$ ASSIGN MAXSHL13.DAT FORO33 
a eee MAXSHL14¢.DAT FQR034 
9 

3! 


RUN PLOTCONV1 


PROGRAM PLOTCONV1] 


THE PURPOSE OF THIS PROGRAM IS TO READ TWO PUNCH=CARD 
FILES FROM PUSLOB AND CONVERT THEM TO A FORMAT WHICH CAN 
BE ACCEPTED BY EASYPLOT. 


DIMENSION VELC100,14),VELO0C100,3),VEL10100,14),LABEL(C14),TAGC14), 
ZALIITEECLS) 
CHARACTERX10 LABEL 
CHARACTER¥40 LTITLE 
COMMON/CSOLVE/DTRECS 
NCURV=14 
DO 10 K=1,NCURV 
READC11,999) LTITLECK),NRECS,DTRECS 
READC11,998) LABELCK), TAGCK) 
READC11,997) CVELCJ,K), J=1,NRECS) 
READC12,999) LTITLECK),NRECS, DTRECS 
READC12,998) LABELCK),TAGCK) 
READC12,997) CVELICJ,K), J=1,NRECS) 
10 CONTINUE 
999 FORMATCAG0,15,1PE11.4) 
998 FORMATCA10,F10.5) 
997 FORMATCIP7E11.4) 
DO 20 L=1,NCURV 
T=0.0 
LL=L+20 
DO 30 M=1,NRECS 
VELQCM,1)=T 
DTRECS1=1000.0XDTRECS 
T=T+DURECS1 
VELOCM,2)=VELCM,L)/412.0 
VELOCM,3)=VEL1(M,L)712.0 
30 CONTINUE 
DQ 40 KK=1,NRECS 
WRITECLL,996) CVELOCKK, JJ),JJ=1,53) 
40 CONTINUE 
WRITECED (99S PETIT leCl) SNRECo EP rRECS 
WRITECLL,994) LABELCL),TAGCL) 
20 CONTINUE 
996 FORMATC3F15.5) 
995 FORMATCAGO,15,1PE11.4) 
994 FORMATCA10,F10.5) 
STOP 
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APPENDIX B 


TYPICAL ELSHOK INPUT CODES--SHELL AND SUBSTRUCTURE 


The input codes which follow are for the resonant 
diaphragm case for both taper and conventional charges. The 
only differences between the two is in the USLOB code where 
the type of charge is specified. All the input codes are 


explained in detail in the ELSHOK users manual [Ref. 4]. 


i. BOSOR4=aNPUT DATA 

The BOSOR4 input code is created by working inter- 
actively with the computer. The prompt following the $ 
symbol appears, and the user supplies the required 
information. The three segments (two endplates and 
cylinder) are modeled separately, and then they are 
"connected" in the global data section. A separate input 
code is required for each circumferential harmonic included 
in the calculation. BOSOR4 output provides the user with 
the in=vacuo free-free modes and matutral frequencies of the 


Shell. 


5 


FILE: BOSSO 


OSORS ATTEMPT oe 
2 


0 
3 
H 
H 
10 
2 
1 
0.0 
0.0 
16.81250 
0.0 
0 
2.500000 
N 


0 
0.0000000E+00 

0 

0 


c 
0.1080000E+08 
0.3200000 
0. 23s 0000E—=05 
0.0000000E+00 

0 

1 


Ztreaa 


45 
3 


1 
16.81250 
0.0000000E+00 
16.81250 
0 ereee 


3 
0.1875000 
N 

6 


2 
3.500000 
53.2500 
103.7500 
154.3750 
205.0000 
254.6250 


batt at feed aed ed 


0.53000000E+03 
7.125000 
7.535500 
2.57 2000 
0.0000000E+00 
“1. 73l250 
0.0000000E+00 
0.7501LE+03 
0.7330000E-03 
0.3000000E+08 
8.250000 
11.69800 
2.750000 
0.0000000E+00 


DATA 


Al 


N=0 


INDIC = ANALYSIS TYPE INDICATOR 


$ NPRT = OUTPUT OPTIONS C1l=MINIMUM, 2=MEDIUM, 3=MAXIMUM) 
ISTRES= OUTPUT CONTROL CO=RESULTANTS, 1=SIGMA, 2=EPSILON) 
NSEG = NUMBER OF SHELL SEGMENTS (LESS THAN 25) 
SEGMENT NUMBER 1 l 1 1 1 1 1 1 
NMESH = NUMBER OF NODE POINTS (5 = MIN.; 98 = MAX.)C€ 1) 
NTYPEH= CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 
NSHAPE= INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
R1 = RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 
el = AXIAL COORDINATE AT BEGINNING OF SEGMENT 
R2 = RADIUS AT END OF SEGMENT. 
fe = AXIAL COORDINATE AT END OF SEGMENT 
IMP = INDICATOR FOR IMPERFECTION CO=NONE, 1=SOME) 
NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
EVAL = DISTANCE (FROM LCEFIMOST SURF. TO UREPERENCE SURE. 


$ 
$ 
$ 
3 
$ 
3 
$ 
$ 
$ 
$ 
$ 
$ 
3 
$ 
$ 
3 
$ 
$ 
3 
$ 
? 
3 
$ 
$ 
$ 
$ 
$ 
9 
3 
$ 
$ 
$ 
> 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
7 
3 
$ 
7 
3 
$ 
$ 
3 
$ 
$ 
$ 
: 
$ 
$ 
3 
3 
$ 
3 
$ 
$ 
> 
> 
> 
$ 
3 
$ 
3 
$ 
$ 


DO YOU WANT TO PRINT OUT RCS), R'CS), ETC. FOR THIS SEGMENT? 
NRINGS= NUMBER (MAX=20) OF DISCRETE RINGS IN THIS SEGMENT 
K=ELASTIC FOUNDATION MODULUS (E.G. LB/INX¥3)IN THIS SEG. 
LINTYP= INDICATOR (0, 1, 2 OR 3) FOR TYPE OF LINE LOADS 
NLTYPE=CONTROL (€0,1,2,3) FOR TYPE OF SURFACE LOADING 
NWALL=INDEX (1, 2, 4, 5, 6+ 7, 8) FOR WALL CONSTRUCTION 

Ee = YOUNG'S MODULUS FOR SKIN 


U = POISSON'S RATIO FOR SKIN 
SM =MASS DENSITY OF SKIN CE.G. ALUM.=.00025 LB-SECHR2/ INK ) 
ALPHA COEFFICIENT OF THERMAL EXPANSION 


oo 


ANRS- CONTROL (0 OR 1) FOR ADDITION OF SMEARED STIFFENERS 
SUR CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

DO YOU WANT TO PRINT OUT THE C(I,J) AT MERIDIONAL STATIONS? 
DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 


SEGMENT NUMBER < 2 c 2 2g 2 ec 

NMESH = NUMBER OF NODE POINTS (5 = MIN.; 98 = MAX.)C 2) 
NTYPEH= CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 
NSHAPE= INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 

RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 


el = AXIAL COORDINATE AT BEGINNING OF SEGMENT 
R2 = RADIUS AT END OF SEGMENT 

ZL = AXIAL COORDINATE AT END OF SEGMENT 

IMP = INDICATOR FOR IMPERFECTION CO=NONE, 1=SOME) 


NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
ZVAL = DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 
DO YOU WANT TO PRINT OUT RCS), R*CS), ETC. FOR THIS SEGMENT? 


NRINGS= NUMBER (MAX=20) OF DISCRETE RINGS IN THIS SEGMENT 
NTYPE = CONTROL FOR IDENTIFICATION QF RING LOCATION (2=Z, 3=R) 
ZCI) = AXIAL COORDINATE OF ITH RING, ZC 1) 

ZCI) = AXIAL COORDINATE OF ITH RING, ZC 2) 

ZCI) = AXIAL COORDINATE OF ITH RING, ZC¢ 3) 

ZCI) = AXIAL COORDINATE OF ITH RING, ZC 4) 

ZCI) = AXIAL COORDINATE OF ITH RING, ZC 5) 

ZCI) = AXIAL COORDINATE OF ITH RING, ZC 6) 

NTYPER= TYPE (0 OR 1 OR 2 OR G OR 5) OF DISCRETE RING NOZt 1) 
NTYPER= TYPE €O OR 1 OR 2 OR @ OR 5) OF DISCRETE RING NO.C 2) 
NTYPER= TYPE (0 OR 1 OR 2 OR @ OR 5) OF DISCRETE RING NO.C 3) 
NTYPER= TYPE (0 OR 1 OR 2 OR G OR 5) OF DISCRETE RING NO.C G) 
NTYPER= TYPE (0 OR 1 OR 2 OR @ OR 5) OF DISCRETE RING NO.C 5) 
NTYPER= TYPE (0 OR 1 OR 2 OR @ OR 5) OF DISCRETE RING NO.C 6) 
E = YOUNG'S MODULUS OF RINGC 1) 

A = CROSS SECTION AREA OF RINGC 1) 

ct = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)¢ 1) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 1) 

LX = PRODUCT OF INERTIAC 1) 

El = RADIAL COMPONENT OF RING ECCENTRICITY (SEE PP.) 700C 1) 
EZ = AXIAL COMPONENT OF RING ECCENTRICITYC 1) 

GJ = TORSIONAL RIGIDITYC 1) 

RM = RING MATERIAL DENSITY CE.G. ALUMINUM=.0002535)¢ 1) 

e = YOUNG'S MODULUS OF RINGC 2) 

A = CROSS SECTION AREA OF RINGC 2) 

1g = MOMENT OF INERTIA ABOUT Y-AXIS CSEE FIG. ON P.70)€ 2) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 2) 

IXY = PRODUCT OF INERTIAC 2) 
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FIDE s0s00 


-2.062500 
0.0000000E+00 
0.9249000E+08 
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0.0000000E+00 
=) 7éi250 
0.0000000E+00 
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eee s 
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Al 
Et = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)€ 2) 
EZ = AXIAL COMPONENT OF RING ECCENTRICITY( 2) 
GJ = TORSIONAL RIGIDITYC 2) 
RM = RING MATERIAL DENSITY CE.G. ALUMINUM=.0002535)¢ 2) 
E = YOUNG'S MODULUS OF RINGC 3) 
A = CROSS SECTION AREA OF RINGC 3) 
LY = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)¢ 3) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 3) 
LEXY = PRODUCT OF INERTIAC 3) 
El = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)¢ 3) 
E2 = AXIAL COMPONENT OF RING ECCENTRICITYC 3) 
GJ = TORSIONAL RIGIDITYC 3) 
RM = RING MATERIAL DENSITY (E.G. ALUMINUM=.0002555)¢ 3) 
E = YOUNG'S MODULUS OF RINGC 4) 
A = CROSS SECTION AREA OF RINGC 4) 
rT = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)¢ 4) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 4) 
IXY = PRODUCT OF INERTIAC 4) 
et = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)¢ 4) 
E2 = AXIAL COMPONENT OF RING ECCENTRICITY( 4) 
GJ = TORSIONAL RIGIDITYC 4) 
RM = RING MATERIAL DENSITY CE.G. ALUMINUM=.0002535)¢ 4) 
E = YOUNG'S MODULUS OF RINGC 5) 
A = CROSS SECTION AREA OF RINGC 5) 
sg = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)¢ 5) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 5) 
IXY = PRODUCT OF INERTIAC 5) 
EX = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)¢ 5) 
EZ = AXIAL COMPONENT OF RING ECCENTRICITYC 5) 
GJ = TORSIONAL RIGIDITYC 5) 
RM = RING MATERIAL DENSITY CE.G. ALUMINUM=.0002535)¢ 5) 
E = YOUNG'S MODULUS OF RINGC 6) 
A = CROSS SECTION AREA OF RINGC 6) 
lY = MOMENT OF INERTIA ABOUT Y-AXIS CSEE FIG. ON P.70)C 6) 
IX = MOMENT OF INERTIA ABOUT X-AXISC 6) 
IXY = PRODUCT OF INERTIAC 6) 
El = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)¢€ 6) 
E2 = AXIAL COMPONENT OF RING ECCENTRICITY( 6) 
GJ = TORSIONAL RIGIDITYC 6) 
RM = RING MATERIAL DENSITY (E.G. ALUMINUM=.0002535)(¢ 6) 


K2ELASTIC FOUNDATION MODULUS CE.G. LB/INX*®3)IN THIS SEG. 
LINTYP= INDICATOR (0, 1, 2 OR 3) FOR TYPE OF LINE LOADS 
NLTYPE=CONTROL (0,1,2,3) FOR TYPE OF SURFACE LOADING 
NWALL=INDEX (1, 2, %, 5, 6, 7, 8) FOR WALL CONSTRUCTION 

E YOUNG'S MODULUS FOR SKIN 

U = POISSON'S RATIO FOR SKIN 

SM =MASS DENSITY OF SKIN CE.G. ALUM.=.00025 LB-SECX¥2/INXXG) 
ALPHA COEFFICIENT OF THERMAL EXPANSION 

ANRS CONTROL (0 OR 1) FOR ADDITION OF SMEARED STIFFENERS 
SUR CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

ARE THERE STRINGERS CPLEASE ANSWER Y OR ND)? 

ARE THERE RINGS CPLEASE ANSWER Y OR ND)? 

K2 =CONTROL (0 OR 1) FOR INTERNAL OR EXTERNAL RINGS 

E2 = RING MODULUS 

U2 = RING POISSON RATIO 

RGMD= RING MASS DENSITY 

IS THE RING CROSS SECTION CONSTANT IN THIS SEGMENT? 

THE RING CROSS SECTION RECTANGULAR CY OR N)? 

ARC LENGTH BETWEEN ADJACENT RINGS CCONSTANT) 
THICKNESS OF RING (CONSTANT) 

HEIGHT OF RING CCONSTANT) 

YOU WANT TO PRINT OUT THE CC(I,J) AT MERIDIONAL STATIONS? 
YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 


SEGMENT NUMBER 3 3 2 3 3 5 3 3 

NMESH = NUMBER OF NODE POINTS (5 = MIN.; 98 = MAX.)C€ 3) 
NTYPEH= CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 
NSHAPE= INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 

RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 

AXIAL COORDINATE AT BEGINNING OF SEGMENT 

RADIUS AT END OF SEGMENT 

AXIAL COORDINATE AT END OF SEGMENT 
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FILE: 60550 
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DATA Al 


DD DADAM P AO OO PMADA DDD AAA DDD AOD DAD AD DAD ADD AMAA ADA DADA AUEAAAUDAUUADUAUUUUDUMOUUH 


IMP = INDICATOR FOR IMPERFECTION (CO=NONE, 1L=SOME) 
NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
ZVAL = DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 
DO YOU WANT TO PRINT OUT RCS), R°CS), ETC. FOR THIS SEGMENT? 
NRINGS= NUMBER (CMAX=20) OF OISCRETE RINGS IN THIS SEGMENT 
KSELASTIC FOUNOATION MOOULUS CE.G. LB/INX¥3)IN THIS SEG. 
LINTYP= INOICATOR (O, 1, 2 OR 3) FOR TYPE OF LINE LOAOS 
NLTYPE=CONTROL (€0,1,2,35) FOR TYPE OF SURFACE LOAOING 
NWALL=INDEX (1, 2, 4, 5, 6, 7» 8) FOR WALL CONSTRUCTION 
E = YOUNG'S MODULUS FOR SKIN 

= POISSON'S RATIO FOR SKIN 
SM =MASS OENSITY OF SKIN (E.G. ALUM.=.00025 LB-SEC#¥2/ IN¥%4) 
ALPHA COEFFICIENT OF THERMAL EXPANSION 
ANRS CONTROL (0 OR 1) FOR AOOITION OF SMEARED STIFFENERS 
SUR CONTROL FOR THICKNESS INPUT (0 OR 1 OR ~-1) 
00 YOU WANT TO PRINT OUT THE CCI,J) AT MERIDIONAL STATIONS? 
DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 


GLOBAL DATA.... 


NOB = STARTING NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS) 
NMINB = MINIMUM NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS) 
NMAXB = MAXIMUM NUMBER OF CIRC. WAVES CBUCKLING ANALYSIS) 
INCRB = INCREMENT IN NUMBER OF CIRC. WAVES (BUCKLING) 

NVEC = NUMBER OF EIGENVALUES FOR EACH WAVE NUMBER 

P = PRESSURE OR SURFACE TRACTION MULTIPLIER 

TEMP = TEMPERATURE RISE MULTIPLIER 

OMEGA = ANGULAR VEL. ABOUT AXIS OF REVOLUTION CRAD/SEC) 


NUMBER OF POLES (PLACES WHERE R=0) IN SEGMENTC 1) 


CONSTRAINT CONDITIONS FOR SEGMENT NO. t 1 i t 
IPOLE = NOOAL POINT NUMBER OF POLE, IPOLEC 1) 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 
IS THIS SEGMENT JOINED TO ANY LOWER=NUMBERED SEGMENTS? 
NUMBER OF POLES (PLACES WHERE R=0) IN SEGMENTC 2) 


CONSTRAINT CONDITIONS FOR SEGMENT NO. a 2 2 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 
IS THIS SEGMENT JOINED TO ANY LOWER~NUMBERED SEGMENTS? 

AT HOW MAY STATIONS IS THIS SEGMENT JOINED TO PREVIOUS SEGS.? 
INODE NODE IN CURRENT SEGMENT CISEG) OF JUNCTION, INOOEC 1) 


JSEG = SEGMENF NO. OF PREVIOUS SEGMENT INVOLVED IN JUNCTION 
JNOOE = NOOE IN PREVIOUS SEGMNT CJSEG) OF JUNCTION 

IUSTAR= AXIAL OISPLACEMENT (CO=NOT SLAVED, L=SLAVEO) 

IVSTAR= CIRCUMFERENTIAL OISPLACEMENT CO=NOT SLAVED, L=SLAVED) 
IWSTAR= RAOIAL DISPLACEMENT CO=NOT SLAVED, 1=SLAVED) 

ICHI = MERIOLONAL ROTATION CO=NOT SLAVED, 1L=SLAVED) 

01 = RAOIAL COMPONENT OF JUNCTURE GAP 

02 = AXIAL COMPONENT OF JUNCTURE GAP 


IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING AND BUCKLING? 
NUMBER OF POLES (CPLACES- WHERE R=0) IN SEGMENTC( 3) : 


CONSTRAINT CONOITIONS FOR SEGMENT NO. 3 3 3 3 
IPOLE = NODAL POINT NUMBER OF POLE, IPOLEC 1) 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 
IS THIS SEGMENT JOINED TO ANY LOWER-NUMBERED SEGMENTS? 

AT HOW MAY STATIONS IS THIS SEGMENT JOINED TO PREVIOUS SEGS.? 


INOOE = NOOE IN CURRENT SEGMENT CISEG) OF JUNCTION, INOOEC 1) 
JSEG = SEGMENT NO. OF PREVIOUS SEGMENT INVOLVEO IN JUNCTION 
JNOOE = NOOE IN PREVIOUS SEGMNT (CJSEG) OF JUNCTION 
IUSTAR= AXIAL DISPLACEMENT CO=NOT SLAVED, L=SLAVED) 


IVSTAR= CIRCUMFERENTIAL OISPLACEMENT (CO=NOT SLAVED, 1=SLAVED) 


IWSTAR= RADIAL DISPLACEMENT CO=NOT SLAVED, I1=SLAVED) 
ICHI = MERIOIONAL ROTATION CO=NOT SLAVED, L=SLAVED) 
D1 = RAOIAL COMPONENT OF JUNCTURE GAP 

02 = AXIAL COMPONENT OF JUNCTURE G 


AP 
IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING ANO BUCKLING? 
DO YOU WANT TO LIST QUTPUT FOR SEGMENTC 1) 
00 YOU WANT TO LIST OUTPUT FOR SEGMENTC 2) 
DO YOU WANT TO LIST OUTPUT FOR SEGMENTC 3) 
00 YOU WANT TO LIST PREBUCKLING RESULTANTS ANO RING FORCES? 


We 


2. ACESNID INPUT DATA 

The ACESNID computer code determines the virtual mass 
array, a quantity producing the late-time contribution of 
the DAA. Only one ACESNID calculation is required for all 
circumferential harmonics of the shell. The bands dividing 
the areas of the shell are identified in the NLEFT, NCYL, 
and NRITE entries, while the fluid properties are specified 


in the RHOFL and VSOUND row. 


FILE: AGEI DATA Al 


VIRTUAL MASS FOR RING-STIFFENED CYLINDER WITH FLAT ENDS, N=0,1 
011410 7 NSTART,NFINIS,NFREQ,NVMASS,NCHECK 

15 51 15 3 0 7 NLEFT,NCYL,NRITE,NWBOSG, NSYMF 

1 6 1000 7 NORDER,NFENDS,NFCENT,NFCMPT,NOMIT 

-5968GE-5 5.833E4 0.005 “ RHOFL,VSOUND, ERR 

Pee Li 7 OUTPUT FLAGS ; 
CPS(1) 

io? Sees, J sopeo, Jn) BEG,J5GEND, JFREND LEFT 

47 / SEFS,JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 

12 “7 SEFS,JSGBEG, JPTBEG, JSGEND, JPTEND, RIGHT 

12 7 SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, LEFT 

47 / SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 


NMWNAQM RRO 
MH eee Ore 
NMNNR NA 
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3. PIFLASH INPUT DATA 

The PIFLASH code takes data from the shell mode files 
from BOSOR4 and the virtual mass file from ACESNID and 
combines them to create a "Shell-fluid file" containing all 
the information concerning the shell and its acoustic medium 
for further calculations. Since only the breathing modes 
are activated in end-on loading, all torsional modes are 
dropped from the shell file. This 1S accomplished in the 


NJUSE rows where he Medes of Table tare Gemenelea. 


PIEGt mee. DATA Al 

2000 0 7 NUMBER,NTORSN,NPTM,NSYMS,NSYMP 

333 3 7 CNNETSGCK),NUSESGCK),K=1,2) 

386.4 7 GRAVITY 

SHELLOO 7 SMF, N=0 

SHELLO1 7 SMF, Ne=l 

19 27 (NJUSEQCIOR =n eho 

123 7 CKORSGCK),K=1,NKORSG) 

294967 9 12 1S 14 19 202i eee coe) 0 ae US epml=0 
1 2/ JUSE, Nel 
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4. SAPIV INPUT DATA 

The substructure is modeled using the standard SAPIV 
code with slight modifications by the developers of ELSHOK. 
The diaphragm 1s modeled using twenty-four plates and 
twenty-five nodes. The location of each node is identified 
in cartesian coordinates. Each plate is 1.06 in thick for 
the resonant case. Nodes eighteen through twenty-five are 
rigidly attached to the shell in the last eight lines of the 
code. The SAPIV code determines the fixed-base modes and 
corresponding natural frequencies as well as the uncon- 


Strained mass and Sti®iness matrices form the substructure. 
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ie Prenuol INPUT DATA 

The PICRUST code takes data from the substructure mode 
file from SAPIV and reorganizes it to facilitate the 
solution of the transient response equations. The location 
of the substructure attachment points to the shell is 


identified for each connecting node. 


FILE: PICRUST1 DATA Al 

bel) 0. io Ome 

Ul lmieeia) 0 1°50 1 7 GUIPUT FLAGS 1-20 
20 0 “ NJUSE,NHWSOB,NHWBAR 

Lome) 0...9 

yee 1) 450 

ZU eae 90.0 

Aimee i155. 0 

2eeceit 130.0 

2c ce. 0 

20 Ne 270.0 

25 2 11 315.0 7 NIPSUB,LBOSEG, LBOSPT,ANGDEG 


0.0 “ DEGROT 
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6. USLOB INPUT DATA 

The two USLOB input codes for the taper and conventional 
charges appear on the following page. USL46 DATA is the 
code for the taper charge. It contains nine discrete data 
points Coon the pressure-time history in seconds and pounds 
per Square inch. There are 1600 time points used in the 
numerical integration of the governing equations with 
information saved every sixteenth point, so only 100 
velocity-time data points are retained for EASYPLOT display. 
The charge is placed along the negative X-axis at 840 
inches. USL47 DATA contains the information for the 
conventional charge. The only difference between the two 
codes 1S a charge identification entry (NCHRG) and the 
charge weight and constants specified for the conventional 


charge instead of the pressure-time history (WCHRG-THEXP). 
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FILE: USL46 DATA Al 


1599 16 29221 / NTIME,NSKIP,NCHRG,NQUAD,NFINE,KOUPLE,NSUBS 

5.0E-05 -840.0 0.0 0.0 7 DELT,XLOAD,RLOAD,SURCUT 

9 / NSHAPE 

0.0 2460.0 0.0002 320.0 0.00053 350.0 0.00243 305.0 0.00294 160.0 

oes 115.0 0.0066 70.0 0.0054 48.0 0.00565 0.0 7 CTCCJ),PCCJ),J=1,NSHAPE) 
0 

0 

4 

i 

1 


ao 


“ DECAY 
1dgdQ0qgq0gd6gd6gd0 QO 

0100000007 QUTPUT FLAGS 1-20 
4 NPTSHL 

ec ll 2 460 12 7 LBOSEG, LBOSP 

4 NPTSUB : 


pie: USLG7 DATA Al 


1599 16 19 221 7 NTIME,NSKIP,NCHRG,NQUAD,NFINE, KOUPLE,NSUBS 
5.0E-05 -840.0 0.0 0.0 “ DELT,XLOAD,RLOAD, SURCUT 
Seems. cootol ters 1.144 35.051351E-5 =.24¢7 7 NCHRG,PZMLT,PZEXP,THMLT, THEXP 


pe 07oc 0000 0 8 

001000000 0 7 OUTPUT FLAGS 1-20 
4 / NPTSHL 

i-imecelle2zece 5 12 / LBOSEG,LBOSPT 
Lo7aNeT SUB 


EOS 


7.  PUSLOB INPUT DATA 

The PUSLOB code produces TEKTRONIX plots of the 
velocity-time histories for specified points on the shell or 
Substructure. It also generates velocity punch card files 


which are converted to a format for EASYPLOT using the 


PLOTCONV1 code. 
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PPee: PUSLS6 DATA Al 


logge to 1 leie7 NIIME NOK IP, NSUBS,NT EK, NGARD 
Seue-U> 1000me).0 7 Delis xMULT > YMULT 
VEPCCI TY PROPILE 
TIME <MSEC) 
VELOCITY CIN/SEC) 
2 NPTSHL-W 
V=enD 1701 
Pe a 00 
V-FWD 3412 
312 -1 1 0.0 7 LBOSEG, LBOSPT,NSPHC, NANG, ANGDEG 
13 NPRTSHAL <U 
V-FWD 1/01 
el Ll lOO: 
voanDs 2702 
eeeeo 1 .0:,0 
V-FWD 2/07 


VeniiDe 2715 
eee 
V-FWD 2720 
22u70 1 
V-FWD 2724 
224 01 
VoRUme7 23 
e26,0.1 
V-FWD 2733 
Zeooe 01 
V-FWD 2737 
eon 
V-FWD 2741 
241 01 
V-FWD 2746 
246 0 1 
VeRNDIS7 12 
Solan. 1 0. 
2 ilione sv 
V-FWD 1701 
ieee =) 10.0 
VSRHD 3712 
312-1 1 0.0 7 LBOSEG,_LBOSPT,NSPHC, NANG, ANGDEG 
1 NPTSUB 
Reeeirr 


20 0 020 0o0 oe 20 92 O& 
> ee o> ee <> Se <> ee => De <> Ee <> es <> ee oD 


/ LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 
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